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, I ABSTRACT 

I ■ Context. The study of symbiotic stars is essential to understand important aspects of stellar evolution in interacting binaries. Their observed 
| population in the Galaxy is however poorly known, and is one to three orders of magnitudes smaller than the predicted population size. 
£C) ■ Aims. IPHAS, the INT Photometric Ho- survey of the Northern Galactic plane, gives us the opportunity to make a systematic, complete search 
£N| ' for symbiotic stars in a magnitude-limited volume, and discover a significant number of new systems. 

Methods. A method of selecting candidate symbiotic stars by combining IPHAS and near-IR (2MASS) colours is presented. It allows us to 
distinguish symbiotic binaries from normal stars and most of the other types of Ha emission line stars in the Galaxy. The only exception are 
£ — ' T Tauri stars, which can however be recognized because of their concentration in star forming regions. 

f^) , Results. Using these selection criteria, we discuss the classification of a list of 4338 IPHAS stars with Her in emission. 1500 to 2000 of them 
!* ' are likely to be Be stars. Among the remaining objects, 1 183 fulfill our photometric constraints to be considered candidate symbiotic stars. The 
. ' spectroscopic confirmation of three of these objects, which are the first new symbiotic stars discovered by IPHAS, proves the potential of the 
. survey and selection method. 

J3; 

Key words. Surveys; Galaxy; stellar content; (Stars:) binaries: symbiotic; Stars: emission-line, Be; Stars: pre-main sequence; (ISM:) planetary 
nebulae: general 
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1. Introduction 

Symbiotic stars are the interacting binaries with the longest or- 
bital periods. They are composed of a compact star, in most 
cases a hot white dwarf, accreting from the wind of a cool 
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giant companion. Part of the giant's wind is ionized by the 
white dwarf, producing the composite spectrum containing 
both absorption features from a cool stellar photosphere and 
emission lines from highly excited ions: these characteris- 
tics originally caused these objects to be named "symbiotic". 
Depending on their near-IR colours, symbiotic stars are di- 
vided into 'stellar' (S) types, if colours are typical of red gi- 
ant branch (RGB) stars, or 'dusty' (D) types, if their near- 
IR emission shows a significant contribution from the warm 
dust known to be typical of evolved asymptotic giant branch 
(AGB) stars. S-types account for around 80% of the sample of 
known symbiotic stars. The presence of an RGB or an AGB 
star also determines the orbital separation at which the sym- 
biotic phenomenon occurs; orbital periods for the S-types are 
between 200 and 6000 days, while they are longer for the D- 
types (>20 yr, but no single robust determination exists so far). 
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For more details on the properties of symbiotic stars, see e.g. 



Corradi, Mikolajewska, & Mahoney (2003) 



A variety of phenomena occur in symbiotic stars that are 
relevant to a number of important astrophysical problems. For 
example, symbiotic stars have been proposed as potential su- 
pernova la progenitors (Munari & Renzini 1992, Hachisu et 
al. 1999). Also, symbiotic stars are excellent laboratories for 
studying (i) thermonuclear outbursts (nova-like accretion in- 
stabilities) under a wide range of conditions (cf. Munari 1997); 
(ii) the powering mechanism of supersoft X-ray sources (cf. 
Jordan et al. 1996); (iii) the collimation of stellar winds and the 
formation of jets (cf. Tomov 2003) and (iv) bipolar (planetary) 
nebulae (cf. Corradi 2003). 

A crucial figure in the discussion of some of these top- 
ics is the total number of symbiotic stars in the Galaxy. 
This is basically unknown. In fact, no systematic search 
for symbiotic stars in the Milky Way has been done so 
far, and the present sample of 173 known Galactic systems 
(Belczynski et al. 2000), plus another 26 suspected ones, is 
mainly the result of occasional discoveries during the study 
of peculiar, variable or erupting stars, or of sparse objec- 
tive prism surveys. This figure should be compared with 
the predicted total number of symbiotic stars in the Galaxy, 
which spans two orders of mag nitude: 3xl0 3 JAllen 19841 , 
3xl0 4 ( |Kenyon et al. 1993|l, 3xl0 5 dMunari & Renzini 1992t , 



and4xl0 3 (Magrini, Corradi & Munari 2003 1. 



IPHAS, the INT Photometric Ha survey of the Northern 
Galactic plane dDrew et al. 2005 ;), gives us the opportunity to 
improve the determination of this basic number. The search for 
symbiotic systems in IPHAS takes advantage of the generally 
strong Ha emission that characterizes this class of objects. In 
this paper, we define our criteria for the detection of symbiotic 
stars using the data from the IPHAS project, complemented by 
near-IR colours from the 2MASS survey. The effectiveness of 
the method is illustrated by the spectroscopic confirmation of 
three new symbiotic stars. 

2. The data 

2.1. IPHAS 

IPHAS is an international collaboration, whose aim is to pro- 
duce a complete, fully photometric, and spatially detailed Ha 
map of the part of the Galactic Plane between latitudes -5° and 
+5° that is visible from the Northern hemisphere. The IPHAS 
observations are obtained using the Wide Field Camera (WFC) 
at the prime focus of the 2.5m Isaac Newton Telescope (INT) 
on La Palma, Spain. The WFC consists of a mosaic of four 
2kx4k EEV CCDs, providing a field of view of 34x34 arcmin 2 
with a sampling of 0".33 per pixel. The IPHAS images are 
taken through three filters: a narrow-band Ha (A c = 6568A; 
FWHM = 95 A) and two broad-band Sloan r and i filters, with 
matched 120, 30, and 10 s exposures, respectively. In this way, 
the magnitude range 13 < r < 20 is covered for point sources 
(the fainter end at 10cr). Pipeline data reduction and data dis- 
tribution are handled by the Cambridge Astronomical Survey 
Unit. The presentation of the survey and further details can be 
found in Drew et al. (2005). 



At the time of writing this article, more than 90% of the 
~1800 square degrees of the northern Galactic Plane to be cov- 
ered by IPHAS has been observed. A first photometric cata- 
logue, containing more than 200 million objects, is about to be 
released (Gonzale z-Solares et al. 2 008). From it, a list of 4853 
Ha emitting stars with r < 19.5 mag has been extracted by 
Witham et al. (2008). This is the sample that we consider in 
this paper for our first search for symbiotic stars within IPHAS. 
Independently, we are employing other methods to select both 
point-like and extended Ha emitters from the IPHAS obser- 
vations (see e.g. Viironen et al. 2008): we will consider these 
additional samples in following papers (note that most symbi- 
otics have a stellar profile, but a small fraction of them might 
be extended in the Ha filter due to their resolved nebulosity, 
see Corradi 2003). 



2.2. Reference samples of known objects 

Unsaturated IPHAS data are available for only four known 
symbiotic stars from the general catalogue of Belczynski et al. 
(2000). They are DQ Ser, V352 Aql, V1413 Aql, and Ap 3- 
1. In order to build a meaningful reference sample of known 
symbiotic stars to be used to define our selection method, addi- 
tional observations were obtained at the INT using the same in- 
strumental setup as for IPHAS. Photometric data of 18 known 
symbiotic stars were taken on June 2, 2004, and April 20, and 
September 17 and 18, 2005. Exposure times were tuned so as 
to avoid saturation of the (generally bright) targets. Reduction 
was done with IRAF. Moreover, existing flux-calibrated spectra 
of 18 symbiotic stars (one of which is in common with the sam- 
ple above), obtained at different epochs, were convolved with 
the WFC filters and instrument response curves to derive Ha, r 
and i magnitudes in the IPHAS photometric system. Where the 
i magnitude cannot be derived from the spectra, we adopted the 
r—i colour from M unari et al. (1992)| All together, we collected 
a sample of 39 known symbiotic stars (29 of the S type and 10 
of the D type) with magnitudes and colours fully consistent 
with those produced by the photometric catalogue of IPHAS. 
This constitutes the primary reference sample to develop our 
selection method. 

We have also investigated other classes of stars and neb- 
ulae which can potentially be confused with symbiotic stars 
according to their IPHAS colours, because of also having Ha 
in emission, or because their spectral type is similar to that of 
the cool component of symbiotic stars. Using IPHAS and ad- 
ditional observations, we have derived Ha, r and i magnitudes 
for 67 planetary nebulae (PNe, Viironen et al. 2008), 79 cat- 
aclysmic variables (Witham et al. 2006), and 518 Mira vari- 
ables. As Be stars are a frequent class of Ha emitting bright 
stars in the Galactic plane (cf. Sect. 14. Il l, they deserve spe- 
cial attention: 18 of these stars in the low-reddening open clus- 
ters NGC 663, NGC 869 and NGC 884, as well as 22 objects 
in the more reddened cluster NGC 7419, turned out to have 
good IPHAS photometric data and were adopted as the com- 
parison samples for this class of stars. Another frequent class 
of Ha emitters in star-forming regions are T Tauri stars; their 
characteristics derived from IPHAS data on Cyg OB2 are be- 
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ing investigated by Vink et al. (2008). Finally, Ha emission is 
also observed in late K to M dwarfs with enhanced chromo- 
spheric and coronal activity (dMe stars). However, their Ha 
emission is generally faint, with line equivalent width^] of a 



few A ( Hawley, Gizis & Reid 1996] ). These are much smaller 
than in symbiotic stars (see Sect. |4 .3b . and for this reason dMe 
stars are not further considered in this paper. 

2.3. Near-IR 2MASS magnitudes 

As symbiotic stars contain a luminous cool giant, near-IR mag- 
nitudes provide information both on their nature and support 
making distinctions between them and other classes of ob- 
jects. The division in the two main groups of S and D types 
was indeed originally made using the J—H and H—K colours 
(Allen & Glass 1974). Recently, the position of symbiotic 
stars in the near-IR colour-colour diagram was discussed by 



Rodriguez-Flores (2006) and Phillips (2007) 



We have therefore extracted from the Two Micron All 
Sky Survey (2MASS) database the J, H and K s mag- 
nitudes of 187 known symbiotic stars, 288 PNe (most 
of them from Ramos-Larios & Phillip s~2005[ ), 95 cata- 
clysmic variables, 1230 Be stars (1148 of which are from 



Zhang, Chen & Yang 2005 1, 487 T Tauri stars (Dahm & Simon 



2005), and 121 Mira variables (Rodriguez-Flores 2006). 

Similarly, we have searched for 2MASS counterparts of the 
4853 Ha emitters in |Withamet al. (2008)| For 4330 of them, 
we found a 2MASS source within 1 arcsec from the IPHAS 
coordinates (note that the two sets of data are calibrated into the 
same astrometric system). Also, the 2MASS counterpart was 
identified for another 8 objects with slightly worse astrometric 
match. This makes a total sample of 4338 Ha emitters with 
2MASS magnitudes. 

It should be noted that the detection limits of IPHAS 
(r<19.5 mag for the objects in Witham et al. 2008) and 
of the 2MASS point-source catalogue {Ks ~15 mag) roughly 
match the characteristic colours of symbiotic stars. We have 
considered a list of 71 known and bright symbiotic stars 
(r<14.5 mag): 70% of them have indeed an observed optical- 
to-near-IR colour (r - Aj)>4.5. On the average, they are pre- 
sumably closer and less reddened than the new symbiotics stars 
that we aim to detect with IPHAS, for which we then expect to 
be able to find a 2MASS counterpart. We are therefore confi- 
dent that adding 2MASS data does not affect the completeness 
of our search for symbiotic stars within IPHAS. 

2.4. Follow-up spectroscopy 

As the present study was progressing, we started a campaign of 
spectroscopic follow-up of the Ha emitters detected by IPHAS. 
Accordingly a dozen candidate symbiotic stars, selected as de- 
scribed in the next sections, were observed at the INT us- 



1 According to the standard definition, the equivalent width of an 
emission line would be negative. However, in order to simplify the dis- 
cussion, in this paper we use a positive sign, so that a larger Ha equiv- 
alent width corresponds to a stronger Ha emission, and to a larger 
r-Hff colour. 



ing the IDS spectrograph, on nights of May 11, June 14, and 
September 9, 2006. Grating R300V was used, which gives a re- 
ciprocal dispersion of 1 .87 A per pixel of the 2kx4k EEV detec- 
tor, and a spectral coverage from 3800 to 8500 A (these figures 
slightly vary from night to night). The slit width was 1.1 arc- 
sec projected on the sky, providing a spectral resolution of 
5 A. Exposure times were 30 min for the two brighter sources 
discussed in Sect. [5] and 2 hours for the faintest and more 
reddened one. Several spectrophotometric standards were ob- 
served during the night for relative flux calibration. Reduction 
was performed using the package onedspec in IRAF. Note that 
the EEV CCD suffers from significant fringing redward of 
~7000 A, which was not possible to remove with the calibra- 
tion frames obtained during those nights. Also, the flux cali- 
bration is somewhat uncertain above 8000 A because of signif- 
icant optical aberrations at the edge of the large format CCD 
used with IDS. 



3. Interpreting the IPHAS and 2MASS 
colour-colour diagrams 

3.1. The IPHAS colour-colour diagram 

Our primary tool for the selection of candidate symbiotic stars 
in the Galactic Plane is the IPHAS r-Ha vs. r—i colour-colour 
diagram, which is presented and thoroughly discussed in Drew 
et al. (2005). It is shown in Fig. [T]for the various classes of 
objects considered. In the diagram, the r-Ha axis mainly in- 
dicates increasing values of the Ha emission line equivalent 
width, while r—i, for normal stars, is essentially a sequence 
of increasing spectral types and/or reddening. In the figure, 
the locus of main-sequence and RGB stars is indicated by the 
solid and dotted lines, respectively (IDrew et al. 2 005). Two se- 
quences are shown, corresponding to reddening values E(B- 
V)=0 (left) and 4 (right), respectively. The reddening vector for 
normal stars, adopted from Howarth (1983) in the same way as 
by Drew et al. (2005), is indicated by the arrow: its length cor- 
responds to 3 mag extinction in V. We note that, as shown by 
Drew et al. (2005), at the higher extinctions probed by IPHAS 
(6 < Ay < 12, roughly) this vector flattens progressively and in 
a mildly spectral-type dependent manner. 

The positions of the different classes of object are indicated 
as follows. Symbiotic stars of the S type are shown as filled cir- 
cles, and the D types as triangles. In the case that several mea- 
surements of the same object obtained at different epochs are 
available, we plot the mean colours as the dispersion of sym- 
biotic stars in the graph is anyway large enough to cover their 
intrinsic variability. PNe are indicated by crosses; cataclysmic 
variables by filled diamonds; Be stars by empty squares; Mira 
variables by dots. In all cases, we show the observed colours. 

As expected, given their emission-line spectra are super- 
posed on relatively weak continua, PNe stand out in the graph 
for their extreme r-Ha colours. Most of them lie close to the 
r-Ha~3.1 limit which is expected for pure Ha emission line 
stars ("ideal" sources for which all the flux in the Ha and r 
bands comes from the Ha line). With smaller Ha excesses than 
PNe (with some overlap), but still above the other classes of Ha 
emitters and normal stars, lie the symbiotic stars. They span a 
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1 2 3 4 5 

r-i 

Fig. 1. IPHAS colour-colour diagram for different classes of objects. The locus of main-sequence and RGB stars is indicated 
by the solid and dotted lines, respectively. Two sequences are shown, corresponding to reddening values E(B-V)=0 (left) and 4 
(right), respectively. The arrow indicates the reddening vector for normal stars: its length corresponds to 3 mag extinction in V. 



significant range in r-i colour. Be stars and cataclysmic vari- 
ables are located closer to the locus of main-sequence stars with 
some moderate Ha excess (except for Nova Aql 1995, which is 
now in the nebular phase and shows a large r-Ha colour). Mira 
variables (which only sometimes show the Ha line in emission) 
are mainly found at the right of the diagram, owing to their 
cool photosphere and frequent large reddening by circumstel- 
lar dust. The class of T Tauri stars is not shown in the figure, 
as they are spread all over the diagram above the sequence of 
normal stars (e.g. Vink et al. 2008). They show a broad range 
of spectral types as well as Ha equivalent widths. The latter 
are in most cases smaller than 500 A (see Fig. 5 in Dahm & 
Simon 2005), which is lower than for the most extreme symbi- 
otic stars. Our strategy to separate them, at least in a statistical 
sense, from the evolved stars on which we are focusing our at- 
tention in this work, will be presented in Sect. 14.41 

The symbiotic stars with the smallest r-Ha colours belong 
to the rarer subgroup of the yellow symbiotics, containing G-K 
giants. Examples are BD-2 13873, AG Dra and LT Del. Also, 
the figure suggests that the D-types have r-Ha colours gener- 
ally larger than the S-types, and more similar to PNe. This was 



already remarked upon in the past (cf. Kenyon 1986), and in 
fact the distinction between PNe and some D-type symbiotic 
stars is tricky (cf. Corradi 2003). However, our present sample 
is limited to ten objects. 

Further information on the location of symbiotic stars in 
the IPHAS diagram, albeit limited to the r-i colour, can be 
obtained from the literature. Mu nari et al. (1992)1 published R 
and I magnitudes for a large sample of symbiotic stars. We have 
transformed these data to the Sloan system for the INT+WFC, 
and in this case we are able to correct for interstellar extinction 
using the data in Whitelock & Munari (1992). The resultant 
range of intrinsic r-i colours spanned by 56 symbiotic stars of 
S type is between -0.1 and +2.8 mag, while the range spanned 
by 19 D-type systems is even larger, from -0.7 to +3.6 mag. 

3.2. The 2MASS colour-colour diagram 

The previous section has shown the potential of the IPHAS 
colour-colour diagram to separate symbiotic stars from the vast 
majority of stars. However, some overlap is present with classes 
of Ha sources whose population in the Galactic Plane is much 
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Fig. 2. 2MASS colour-colour diagram for different classes of objects. Symbols are the same as in Fig. [1] albeit smaller owing to 
the larger number of objects. In addition, T Tauri stars are indicated by empty circles. Like in Fig.Q] the locus of main-sequence 
and RGB stars is indicated by the solid and dotted lines, respectively, for reddening values E(B-V)=0 (lower-left sequence) and 
4 (upper-right sequence). The arrow shows the reddening vector for normal stars corresponding to 3 mag extinction in V. 



larger than that of symbiotic stars. In addition, mixing between 
classes is inevitably raised in the presence of photometric er- 
rors. 



The 2MASS J-H vs. H-K$ colour-colour diagram pro- 
vides an additional resource to refine the selection of symbi- 
otic stars (Rodnguez-Flores 2006). The various classes of ob- 
jects, as well as normal main-sequence and giant stars (Bessel 
& Brett 1988), are indicated in Fig. [2] with the same symbols 
as in Fig. Q] In addition, we plot here as empty circles 104 
classical T Tauri stars from Dahm & Simo n (2005)| with an Ha 
equivalent width larger than 40 A (roughly corresponding to 
the smallest Ha equivalent width shown by the known sym- 
biotic stars in Sect. 13.11 ). All these T Tauri stars belong to the 
young cluster NGC 2264, which has a low foreground extinc- 
tion (Ay = 0.22 mag), and is rich in pre-main sequence Ha 
emitters (Dahm & Simon 2005). For all objects we plot the ob- 
served colours. As for the IPHAS diagram, the reddening vec- 
tor is computed from Howarth (1983) for the 2MASS filters J, 
H, and K s : the adopted shifts for a reddening of 1 mag in V are 
0.1 13 and 0.069 for the J-H and H-K$ colours, respectively. 



The bulk of PNe occupy a distinct locus at the bottom of 
the near-IR diagram (see also Ramos-Larios & Phillips 2005), 
to the right of normal stars. Most of the cataclysmic variables 
and Be stars are instead clumped closer to the locus of main- 
sequence stars. In particular, note the well-defined sequence at 
the bottom of the diagram displayed by Be stars, which extends 
to the right side of the unreddened main-sequence stars. We 
will come back to this important feature in Sect. 14.11 



S-type symbiotic stars, according to their original defini- 
tion (Allen & Glass 1974), have colours typical of RGB stars, 
with an upper-right tail due to reddening. They form a compact 
cluster in the diagram, except for few objects. The most ex- 
treme case is UV Aur (a carbon star), which is displaced to the 
right of the S-types sequence, toward the broad region which 
is occupied by T Tauri stars, D-type symbiotic stars, and some 
extreme Mira variables. The location of D-type symbiotics is 



modeled by Phillips (2007) with a variable combination of stel- 
lar and warm dust continuum emission, plus some extinction. 
Three objects, AS 201, V471 Per, and StHA 190, have signifi- 
cantly smaller J-H colours than the bulk of symbiotics. Again, 
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Fig. 3. The 4338 Ha emitters by |Witham et al. (2008)| in the 
IPHAS (top) and 2MASS (bottom) colour-colour diagrams. 
The loci of main-sequence and giant stars and the reddening 
vectors are also indicated (see Fig. [TJ for a detailed explana- 
tion). 



these are yellow symbiotics, which tend to occupy the bottom 
part of the locus of symbiotic stars, as in the IPHAS diagram. 

Some PNe fall in the region of D-type symbi- 
otic stars. These are examples of possible misclassi- 
fied PNe, which might instead be symbiotic stars with 
extended nebulae as frequently suggested in the lit- 
erature ( ICorradi 19951 



Schmeja & Kimeswenger 2001 



Santander-Garcfa, Corradi & Mampaso 2007|l. 



Finally, T Tauri stars form a broad sequence running par- 
allel to, and on right side of, S-type symbiotics, that extends 
well into the locus of the D-type symbiotic stars. Note that 
T Tauri stars with smaller equivalent widths than those dis- 
played, including the so-called weak-line T Tauri stars, would 
be clumped on the left-bottom end of this sequence (Dahm & 
Simon 2005). The significant overlap with symbiotic stars in 



both the 2MASS and IPHAS diagrams, makes T Tauri stars 
the most frequent "contaminants" in our search for symbiotics 
stars. 

4. Application to the list of IPHAS Ha emitters by 



Witham et al. (2008) 



The analysis of the IPHAS and 2MASS colour-colour diagrams 
presented in the previous section allows us to discuss the nature 
of the 4338 IPHAS Ha emitters by |Witham et al. (2008)| with 
2MASS counterparts (Sections 12. 1 1 and [ 2.3b . Their location in 
the colour-colour diagrams is shown in Fig. [3] Looking at the 
2MASS diagram, there are some striking similarities with the 
plot of known objects in Fig. [2] First, the 2MASS diagram of 
Fig. |3] shows a main concentration of sources at its bottom, to 
the right of the unreddened main-sequence stars. They form a 
dense band with a slightly lower inclination than our adopted 
reddening vector. This is obviously similar to the sequence de- 
fined by known Be stars in Fig. [2] Second, another inclined 
sequence of objects can be seen in Fig. [3] at higher J-H val- 
ues, in the region where T Tauri stars (frequent in the Galactic 
Plane) and symbiotic stars (presumably less numerous) should 
be located. 

We therefore start our discussion by considering what ap- 
pears to be one of the most represented classes of objects in the 
list of |Witham et al. (2008")] the Be stars. 

4.1. Be stars 

In order to further investigate the similarities between Fig. [2] 
and|3] we have studied in more detail the location of Be stars in 
the colour-colour diagrams. Zhang, Chen & Yang (2005) pro- 
vide the reddening of 1 148 individual objects: this allow us to 
identify the locus of unreddened Be stars in the 2MASS di- 
agram, and to define a corresponding "selection" box. These 
are indicated by the dots and the solid line, respectively, in the 
lower-left panel of Fig. [4] When shifted to the reddening of the 
clusters NGC 663, NGC 869 and NGC 884 (A v ~2 mag, dotted 
line) and NGC 7419 {A v ~6 mag, dashed line), such a box in- 
cludes the corresponding data points for these clusters (empty 
and filled squares, respectively). 

Similar consistent results are obtained for the IPHAS 
colour-colour diagram. There, if we define the locus of Be stars 
from the position of the objects in the four open clusters con- 
sidered, and correct for the corresponding reddening, we can 
determine a zero-reddening selection box for this class of ob- 
jects, which is shown as a solid line in the upper-left panel of 

Fig.a 

We have then done the exercise of extracting can- 
didate Be stars from the 4338 IPHAS Ha emitters in 
With am et al. (2008)| An object is considered to be a good can- 
didate Be star if it falls inside the selection boxes in both the 
IPHAS and the 2MASS diagrams. The combination of the se- 
lection boxes in the two diagrams, shifted so as to consider 
a range of reddening values from Ay-Q to Ay=8 mag (repre- 
sented by the sequence of boxes in the right panels of Fig. HJ, 
results in a total of 2035 candidate Be stars (also plotted in 
Fig. |4}. The majority of them would have reddening values 
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Fig. 4. Left: location of known Be 
stars in the IPHAS and 2MASS 
diagrams. The dots mark dered- 
dened colours for the 1 148 Be stars 
from Zhang, Chen & Yang (2005). 
Open squares are the reddened data 
for the open clusters NGC 663, 
NGC 869 and NGC 884, and filled 
squares similar data for NGC 7419. 
The solid line is the adopted zero- 
reddening selection box for Be 
stars. Right: the 2035 candidate Be 
stars in the list of 4338 Ha emit- 
ters by |Witham et al. (2008)| The 
alternate solid/dotted boxes show 
the selection boxes (see text) for 
A, ,=0,2,4,6,8 mag, from bottom-left 
to top-right, respectively. 



between 3 and 5 magnitudes in V. Objects at zero redden- 
ing are virtually absent, presumably because such objects are 
in the main nearby and thus bright and saturated in IPHAS. 
Only for significant reddening values (i.e. at the largest J—H 
colours, cf. the lower-right panel of Fig. |4] with Fig.[2]i, mixing 
of the Be stars candidates with low-reddening T Tauri stars is 
possible. An idea of the magnitude of the contamination can 
be gained considering that only 4% of the 450 T Tauri stars 
in |Dahm & Simon (2005) including those with the lowest Ha 
equivalent widths, have J-H< 0.55. Adopting this figure, and 
the worst possible (and unlikely) scenario that all sources with 
J—H> 0.55 in our complete list of 4338 objects are T Tauri 
stars, it would be concluded that only one hundred objects of 
this class would be expected with a lower J—H colour, at vari- 
ance with the more than 1600 Be star candidates (80% of our 
list) lying below this colour limit. 

A smaller amount of mixing is expected with other, less 
frequent classes of objects, like CVs, compact PNe, and sym- 
biotic stars. We conclude that somewhere between 1500 and 
2000 Ha emitters in \Witham et a l. (2008) are likely to be Be 
stars, which would therefore be the largest class of objects in 
this list of IPHAS Ha stars. 

4.2. Planetary nebulae 

In Sect. 13.11 we remarked on the overlap in the IPHAS di- 
agram between PNe and symbiotic stars, especially those of 
the D-type. On the other hand, the 2MASS colours are gen- 
erally different, allowing a neat separation of the two classes 
of objects (with notable exceptions which might be related to 
misclassification). A search for compact PNe from the list of 



|Witham et al. (2008)1 resulting in a few new candidates, and a 
confirmation of their nature by spectroscopic observations, is 
presented elsewhere (Viironen et al. 2008). 

4.3. Symbiotic stars 

In a similar way as done for Be stars, based on the properties of 
our reference samples, we define selection boxes in the colour- 
colour diagrams to be used for our search for symbiotic stars in 
the Milky Way. 

In the IPHAS diagram, the limited number of symbiotic 
stars of the D-type in the reference sample, as well as their 
significant overlap with the S-types (with possible differences, 
however, as outlined in Sect. 13. U . prevent us from defining sep- 
arate boxes for the two classes. We prefer instead to use the 
same selection criterion, leaving to the near-IR colours the task 
of separating the two classes according to the original defini- 
tion. Our selection box for symbiotic stars (Fig. [5] upper panel) 
is then defined as follows. In the r—i axis, we allow the range 
of intrinsic colours spanned by known objects as discussed at 
the end of Sect. 13.11 and extend it to the right side of the graph 
to allow for the significant extinction that is expected when ob- 
serving through the Galactic Plane. 

In the vertical direction, we define a lower limit which is 
an inclined line, parallel to the reddening vector (Sect. 13. U . 
and defined by the formula (r - Ha) > 0.25 ■ (r - i) + 0.65, 
which includes all known symbiotic stars in the samples of 
Sect. 13.11 except for the yellow symbiotic BD-213873. This 
limit roughly corresponds to Ha equivalent widths of 50 A 
(Drew et al. 2005). We decided to fix the lowest side of the 
selection box at this limit in order to avoid significant mixing 
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Fig. 5. Top: the selection box for symbiotic stars in the IPHAS 
colour-colour diagram, superimposed on the reference sample 
and the locus of main-sequence and RGB stars (see Fig. Q] for 
details). Bottom: the selection boxes for the S-type (continuous 
line) and the D-type (dashed line) symbiotics in the 2MASS 
diagram. Symbols as in Fig. [2] 



with Be stars (see Sect. 14. 11 1. The upper limit of the selection 
box for symbiotic stars is the r-Ha value for pure Ha emitters, 
namely r-Ha~3 . 1 . 

In the 2MASS diagram, the large number of objects avail- 
able allows us to define separate boxes for the S and D-type 
symbiotic stars. For the S-types, we first set the lower limit for 
the J-H colour which includes all known objects. Then the left 
and right limits for the H-Ks colours were chosen to run paral- 
lel to the reddening vector. They are chosen to include the vast 
majority of S-type systems, except for nine objects (7% of the 
whole sample) that are detached and on the right side - these 
will, in any case, be included in the selection box for D-types. 
We have decided to keep the right limit for S-types (the most 
frequent class of symbiotics) as leftward as possible, in order 



to minimize mixing with T Tauri stars (even if at the expense 
of the D-types). The selection box for S-type symbiotic stars is 
indicated by the solid line in the lower panel of Fig. [5] For the 
D-types, we use the same procedure, with the further assump- 
tion of setting the same lower limit for the J-H colour as for 
the S-types, and the H-K$ left limit so as to have contiguous 
boxes with no gap or overlap. The selection box for the D-type 
symbiotics is indicated by a dashed line in Fig. [5] 

In this way, the three yellow symbiotics AS 201, StHA 190, 
and V471 Per, fall out of both boxes in the 2MASS diagram. 
But they are located in the very populated region of the diagram 
where not only main sequence stars are found, but also other 
Ha emitters like Be stars and cataclysmic variables. With so 
much mixing, a search for this kind of object becomes very 
difficult. 

When our selection criteria are applied to the list of 
IWitham et al. (2008)| 337 sources fall in both the IPHAS and 
2MASS selection boxes for the S-types (Fig. left), and 846 
fulfill the criteria to be considered D-type candidates (Fig. [6] 
right). The complete list of all these candidates is reported in 
Tab.Q] which is only available electronically. In the 2MASS di- 
agram, a good number of objects fall in the lower right part of 
the selection boxes for S-type symbiotic stars, while most of the 
objects in the D-type box fall in its left-bottom side. These are 
the regions where T Tauri stars are also expected to be found, 
indicating that significant mixing of the two classes of objects 
is likely to be present in the candidate list (Tab.[T]i. 

Our selection method recovers the three known sym- 
biotic stars included in the list of |Witham et al. (2008)] as 
well as two suspected ones. Out of the 29 objects listed by 
|Witham et al. (2008)1 as known young stellar objects, 13 are 
included in the list of D-type candidate symbiotics, and 3 
in the list of S-type candidates (note that none of them is 
classified as a Be star). This confirms that the mixing with 
young stellar objects is more severe for the D-types. Among 
the 9 known or newly discovered PNe included in the list of 
|Witham et al. (2008)] only one enters the list of D-type candi- 
dates. Its nature will be discussed in Viironen et al. (2008). 

Note that if the lowest edge of the IPHAS selection box is 
lowered so as to include all known symbiotic stars in our sam- 
ples, then another 152 candidates of S-type, plus 330 of D-type, 
would be added. At the same time, the confusion with other 
classes of object would also increase. In this respect, we stress 
that our aim is not to define absolute limits for the colours of 
symbiotic stars, but only to provide a practical way of select- 
ing new candidates from the IPHAS survey, guided by the best 
reference samples that we were able to build from the litera- 
ture, IPHAS itself, and from additional observations obtained 
in backup time during our observational campaign. As demon- 
strated below, the proposed selection method seems indeed to 
be promising. 

4.4. Clustering as an additional criterion to separate 
young stars from symbiotics 

Observations targeting a narrow band of the Galactic Plane 
must include a large number of star-forming regions, young 
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Fig. 6. Left: The 337 objects with 
IPHAS and 2MASS colours fulfill- 
ing our selection criteria for S-type 
symbiotic stars. The locations of 
the three new symbiotic stars dis- 
covered by IPHAS, for which we 
present spectroscopic confirmation 
in Sect. [5] are indicated by the cir- 
cles. Right: the 846 candidate D- 
type symbiotic stars. 



clusters and associations. There, numerous objects, especially 
Be and T Tauri stars, are Ha emitters. We have shown that 
while Be stars can be separated from symbiotic stars because 
of their smaller r-Ha and J—H colours, T Tauri stars over- 
lap with symbiotics in both the IPHAS and 2MASS diagrams. 
Therefore, they must be regarded as the most serious "con- 
taminants" in our search for symbiotic stars using IPHAS and 
2MASS data. 

Our strategy for tackling the problem is to consider the spa- 
tial distribution of the Ho- emitters selected by IPHAS. T Tauri 
stars are expected mainly to be concentrated in young clusters. 
On the other hand, symbiotic stars, which belong to an older 
Galactic stellar population (bulge/thick-disc, Munari & Renzini 
1992), should be more isolated Ha sources (except for projec- 
tion effects through the Galactic Plane). 

A first step in exploring this issue was taken by mea- 
suring the degree of clustering of the objects in the list of 
IWitham et al. (2008)[ For each of the 4853 Ha emitters, the 
mean distance to the n most nearby objects (with n from 1 to 
12) was computed. We show in Fig. [7] the distribution of the 
distances for n = 4, that we find to be a good compromise be- 
tween having a sufficient number of neighbours to detect the 
existence of a group, while avoiding the limited statistics that 
the requirement of a large number of neighbours might suffer 
from, given the size of the global sample. 

The 29 objects listed by With am et al. (2008)| as young stel- 
lar objects, as well as 20 known, suspected or confirmed new 
PNe (Viironen et al. 2008) and symbiotic stars (to be presented 
in paper II) were then considered. They confirm the expected 
trend: young stars (the vertical short lines at the top of Fig. |7]i 



tend to occur with smaller angular separations than do PNe and 
symbiotic stars (the filled circles). Therefore, this additional 
parameter seems to be a useful one to distinguish, at least in 
a statistical sense, young stellar objects from stars belonging 
to older Galactic populations like symbiotic binaries. We ex- 
pect that the significance of this parameter will improve fur- 
ther as the sample of Ha emitters detected by IPHAS grows in 
number; this would allow a better definition of clustering in the 




20 40 60 80 100 

mean distance from 4 neighbours [arcmin] 

Fig. 7. Histogram of the mean distance to the four most 
nearby objects for each of the 4853 Ha emitters by 
|Witham et al. (2008)1 29 young stellar objects are indicated by 
the short lines at the top of the diagram, while 20 PNe and sym- 
biotic stars are indicated by the dots. 
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area observed by IPHAS. For the time being, we add the mean 
distance to the four nearest neighbours as a further datum in 
the list of candidate symbiotic stars in Tab. [TJ the weight to 
be given to this parameter is left to the discretion of the user. In 
any case, we note that about 50% of both the S-type and D-type 
candidates in Tab.[T]have a mean distance to their four nearest 
neighbours that exceeds 15 arcmin: these candidates are our 
preferred initial targets for spectroscopic follow-up. 

5. Spectroscopic confirmation: the first symbiotic 
stars discovered by IPHAS 

As an illustration of the potential of IPHAS and the selection 
method proposed in the previous sections, we present here the 
spectra of the first three new symbiotic stars discovered by 
IPHAS. Data for more stars and a deeper analysis of the in- 
dividual objects will be presented in paper II. 

In the following, we use both the standard nomenclature for 
IPHAS point sources ("IPHAS J" followed by the J2000 coor- 
dinates, see Drew et al. 2005), as well as the name "IPHAS-Sy 
nn" to indicate the ««th symbiotic star discovered by the sur- 
vey. 

The INT spectra are shown in Fig. [8] and the location of 
the objects in the IPHAS and 2MASS diagrams is indicated 
by open circles in Fig. [6] The three objects have large r-Ha 
colours and fall in the selection box for S-type symbiotics. In 
the latter, two of them have colours characteristic of strongly 
reddened red giants, as confirmed by the INT spectra. 

5.1. IPHAS J1 85704.44+002631. 7 (IPHAS-Sy 1) 

The IPHAS and 2MASS magnitudes of IPHAS-Sy 1 are: 
r=18.28, /=16.03, Ha=15.59, 7=10.42, //=8.86, and 7^=8.09. 
No objects are listed in the SIMBAD database, within a ra- 
dius of 2 arcmin about the source coordinates, and the nearest 
known radio HII detection is over 30 arcmin away. There is no 
evidence to link this object to a known star-forming region. 

The Ha equivalent width is 800 A, which corresponds to 
r-Ha=2.69. This is the faintest objects of the three in the opti- 
cal (but not in the near-IR), and the one with the highest 7-77 
colour. This can be mainly ascribed to reddening, as shown by 
the steep decrease of the flux below 5000 A (no emission is de- 
tected blueward of the [OIII]5007 A line) and the large r - K s 
colour amounting to 10.2 magnitudes. The symbiotic nature of 
IPHAS-Sy 1 is indicated by the simultaneous presence of high 
excitation emission lines of [OIII] and [Arlll], and the TiO ab- 
sorption bands of a M-type star (Fig. [8] upper panel). 

5.2. IPHAS J1 90924.64-01 091 0.2 (IPHAS-Sy 2) 

IPHAS-Sy 2 has the following magnitudes: r=17.16, /=15.26, 
Ha=14.79, 7=11.44, 77=10.29, and 7^=9.80. Its Ha equiva- 
lent width is 650 A (r-Ha=2.38). 

In SIMBAD, it is indicated as a possible PN (K 4-17), but 
the spectrum described by |Stenholm & Acker (1987j| does not 
allow its real nature as a symbiotic star to be identified. This 
is instead revealed by our INT spectroscopy, which captures a 
rich set of high excitation emission lines ([Nelll], [OIII], Hell, 



and even [FeVII]), accompanied by the red continuum of an 
M4 giant (Fig. [8] middle panel). No star formation regions are 
catalogued in the vicinity of IPHAS-Sy 2. 

5.3. IPHAS J20251 0.58+435233.0 (IPHAS-Sy 3) 

Its magnitudes are: r=15.34, /=13.34, Ha=13.42, 7=9.94, 
77=8.45, and 7Cs=7.73. The Ha equivalent width is 250 A 
(r-Ha=1.92). 

IPHAS-Sy 3 is listed in SIMBAD as a misclassi- 
fied planetary nebula (PN K 3-59), but the spectrum of 
Sabbadin , Falomo & Ortolani (1987)1 on ly reveals weak Ha 
line emission superimposed on a very red continuum. Our ob- 
servations show the composite spectrum defining a symbiotic 
star, with the simultaneous presence of the continuum of an 
early M star and high excitation emission lines. 

In projection, the star is located in the Cygnus X com- 
plex, a little over 3 arcmin from the HII region DWB 137 
(Dickel, Wend ker & Bieritz 1 969). The vicinity of young stars 
might cast some doubt on the nature of IPHAS-Sy 3 as a sym- 
biotic star. Even if it shows some of the typical spectral sig- 
natures of these interacting binaries, such as the TiO bands to- 
gether with highly ionized Hell and [OIII] in emission, some 
extreme examples of T Tauri stars also present them (see e.g. 
Beristain, Edwards & K wan 200l] l. 

However, there are several arguments pointing toward 
IPHAS-Sy 3 as most likely being a symbiotic star. First, its 
spectrum lacks evidence of other features typical of T Tauri 
stars: e.g. the strongest resonance Li doublet at 6708 A in ab- 
sorption; any of Fell or Nal lines in emission; apparent broad- 
ening of the Hel emission line profiles. Second, T Tauri stars, 
being closer to the main sequence (MS), have radii just 2- 
3 times larger than their MS descendants, whereas symbiotic 
stars must contain a red giant with radii two orders of magni- 
tudes larger than MS values. Playing devil's advocate, we can 
see if a typical T Tauri star radius works for IPHAS Sy 3 if 
placed at a plausible distance to be in a young Cygnus-X clus- 
ter. The distance to the young cluster containing the Herbig Be 
star VI 685 Cyg, some 3 degrees away on the sky, is given as 
980 pc (Davi eTet al. 20011) . with the other nearby HII regions 



in Cygnus X being mostly at 1-2 kpc ( Straizys et al. 1993 1. For 
a distance of ~1 kpc, an M0 spectral type (compatible with 
the strength of the TiO and VO absorption bands), a redden- 
ing E(B-V)~2-2.5 (estimated from HI and Hel emission line 
ratios, the J-Ks colour, and a fit of the optical and near-IR 
spectral energy distribution), we apply the Barnes-Evans rela- 
tion (e.g. Beuermann, Baraffe, & Hauschildt 1999) to derive a 
radius ~30 R G . This is too large for a T Tauri star, but consistent 
with the presence of a red giant. The estimated radius increases 
if the distance is larger than 1 kpc. 

6. Summary and perspectives 

Searching for a relatively rare and old family of stars in a nar- 
row band in the Galactic plane, where young stellar objects 
dominate the population of Ha emitters, might seem a hard 
challenge. Our motivation to embark in such an enterprise is 
the poor knowledge of the total population of symbiotic stars 
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Fig. 8. Spectra of the first three new symbiotic stars discovered by IPHAS. The region around 7600 A with the strong oxygen 
atmospheric absorption band is not plotted. The main low and high excitation emission lines and the heads of the TiO bands, 
which demonstrate the symbiotic nature of the objects, are labeled. The oscillations in the very red part of the spectra are due to 
fringing in the detector. 



in the Galaxy: the IPHAS survey gives us the opportunity to 
tackle this problem by performing, for the first time, a com- 
plete search of these objects in a magnitude limited volume. 

As symbiotic systems contain a cool giant star, the near-IR 
data from the 2MASS survey have been added to the IPHAS 
photometry. Then, a discussion of the location of the different 
classes of Ha emitters in the combined IPHAS and 2MASS 
colour-colour diagrams has been presented. This allows us to 



define selection criteria for symbiotic stars which separate them 
from the vast majority of normal and Ha emitting stars. 

The only exception are T Tauri stars. They overlap with 
symbiotic stars in the IPHAS diagram. In the 2MASS one, 
they form a sequence very close to (and partially overlap- 
ping with) the main subclass of symbiotic stars (S-types), and 
fully overlap with the other, less frequent group of symbiotics 
(the D-types). In fact, even when a spectrum is available (see 
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Sect. 15. 3K it is sometimes difficult to distinguish a symbiotic 
star from those T Tauri stars which display high excitation 
emission lines. As an additional criterion to separate symbi- 
otic from T Tauri candidates in our survey, we have considered 
the spatial distribution of the Ha emitters detected by IPHAS. 
T Tauri stars are expected to be found in groups in star forming 
regions. In contrast, symbiotic stars should be more isolated 
objects. With this aim, we have introduced a "clustering" pa- 
rameter which is the mean distance to the four most nearby Ho- 
emitters in the IPHAS list that we have considered (Witham et 
al. 2008). 

We expect that all these criteria for the search for symbiotic 
stars within IPHAS will improve as soon as additional infor- 
mation is gathered from the survey itself. A better fix on the 
number of Ha emission line stars above specified, more exact- 
ing equivalent width thresholds, might improve our estimates 
of the clustering parameter - or further follow-up spectroscopy 
of the selected candidates can improve our understanding of the 
IPHAS and near-infrared colour-colour planes. 

With this second aim in mind, we have started an intensive 
spectroscopic campaign. With it, we hope that the number of 
symbiotic stars in the IPHAS area, so far limited to only 1 1 
objects listed in the catalogue of Belczynski et al. (2000), will 
significantly increase. The confirmation in this paper of three 
new systems out of a small sample of candidates observed in 
the first nights of the spectroscopic follow-up, seems to be a 
good start in this direction. 
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Table 1. 1 183 candidate symbiotic stars extracted from the list 
of |Witham et al. (2 008) using the photometric constraints de- 
fined in Sect. 14. 31 Columns contain: the coordinates, the IPHAS 
and 2MASS magnitudes and their errors, the mean distance d\ 
(in arcmin) from four neighbours (see Sect. 14. A\ . Tab. la con- 
tains the candidate S-type symbiotics, and Tab. lb the candidate 
D-types. The Table is available only in electronic form via the 
CDS. 



